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Abstract 
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1 Introduction 



There are several methods to construct the Darboux matrix (which gener- 
ates soliton solutions) ^HHI3I3JIII1O30IE])- However, these methods are 
technically difficult when applied to the matrix versions of the spectral prob- 
lems which are naturally represented in Clifford algebras [HI I1UI I12| . Some 
of these problems are avoided in our recent paper ^3] . In the present paper 
we develop the ideas of ^3] in the matrix case. We extend our approach on 
the multisoliton case and consider the reduction problem and the discrete 
case. We also show that our approach, although different, is to some extent 
equivalent to the standard dressing method. We compare our method with 
the Zakharov-Shabat approach ^ and the Neugebauer-Meinel approach 

H3Q2]. 

We consider the spectral problem 

*„ = £V*, = l,...,m) (1) 

(with no assumptions on J7„ except rational dependence on A) and the Dar- 
boux transformation 

# = , (2) 

which means that 

= , (3) 

where and Un have the same rational dependence on A (U^ and ^ are 
n x n matrices but our approach works well also in the Clifford numbers 
case [T3]^. 

The construction of the Darboux transformation is well known (espe- 
cially in the matrix case) [7lll4j. The first step is the equation for D resulting 
from ©,© and ©: 

D,^+DU^ = U fl D . (4) 



In our erlier paper |T3] we proposed the following procedure. We assume 
that there exist two different values of A, say A+ and A_, satisfying 

D 2 (\±) = . (5) 

Denoting ^(A±) = D(\±) = D±, evaluating @ at A = A± and multi- 
plying @ by D± from the right, we get: 

D±, fl D ± + D±U fl (\±)D± = . (6) 
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We assume that ^(A±) are invertible (which is obviously true in the generic 
case). It is not difficult to check that D± given by 

At = <p±V±d±^ ,4 = 0, (7) 

(where d± = const and <p± are scalar functions) satisfy equations (J5J), ©. 
Assuming that D is linear in A, i.e., 

D(\)=A + A 1 X, (8) 

we can easily express Aq , A% by D± to get 

d(x) = + j^r^-d-*- 1 ■ ( 9 ) 



2 One-soliton case and the Zakharov-Shabat ap- 
proach 

We confine ourselves to the case linear in A (see ©). The condition (JSJ) can 
be easily realized if 

D 2 (X) =ct(A-A+)(A- A_)J (10) 

where u ^ is a constant, A+ ^ A_ and / is the identity matrix. The 
identity matrix will be sometimes omitted (i.e., for a 6 C we write al = a). 
In the case ([10[) from © and Q it follows that 

D + D_ + D_D + = -a{\ + -\-) 2 . (11) 



Lemma 1 D of the form {SJ) satisfies hl(J\) if and only if n is even and 

D=N(\-\+ + (\ + -\-)P) (12) 

where the matrices AT and P satisfy 

P 2 = p , TV 2 = a , MPN- 1 = I-P . (13) 

In this case the Darboux matrices and are equivalent. 

Proof: We denote TV :— A\. From JSJ) we get 
D 2 {\) = Al + (A a Af + J\fA )\ + TV 2 A 2 , 
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i.e., D 2 (X) is a quadratic polynomial. It is proportional to the identity matrix / 
(compare JTUJ)) iff 

J\f 2 = a , A JV + JVA = -a(X+ + A_) , A 2 = <tA + A_ . (14) 
Multiplying the second equation by AfAo we get 

a 2 A+A__ + (NA Q ) 2 + a(X+ + \_)J\fA Q = . 
Hence (AfAo + <r\ + )(AfAo + <tA_) = 0, and, denoting Q := MAq + <tA + , we have 

Q 2 = {X+-X-)aQ 

which means that Q — (A + — A_)erP, where P 2 = P. Therefore, taking into account 
TV 2 = (J, we get ()12|l. Now, we take into account the third equation of (|14J) . First, 
A\P = aX+X-P yields A_(A+ - X-)MPMP = 0. Then the equation A 2 = aX+X- 
is equivalent to A+(A+ - A_)(cr(7 - P) - MPN) = 0. Therefore AfPAf- 1 = I -P. 
This equality means that kerP = AA _1 ImP which implies dimkerP = dimlmP. 
Thus n is even which complets the proof. □ 

The case A + = A_ can be treated in a similar way and it leads to the 
nilpotent case 0: 

D = M(X -X+ + M) , M 2 = , AT 2 = a , M = -NMN' 1 . 

Our method is closely related to the standard dressing transformation 
^ El E| • The Darboux matrix can be rewritten as 

D = {X -\ + )n(i + ^^-P^ . (15) 

We recognize the standard one-soliton Darboux matrix in the Zakharov- 
Shabat form [7JE1- We point out that usually one considers the Darboux 
matrix T> = (A — A+) _1 P> which is equivalent to D given by (|12[) because the 
multiplication of D by a constant factor leaves the equation Q invariant 
|16j . M is known as the normalization matrix and P is a projector expressed 
by the background wave function: 

kerP = ^{X+)V ker , imP = ^{X-)V im , (16) 

Vk er and Vim are some constant vector spaces, A + and A_ are constant com- 
plex parameters. The last constraint of (jl<-i|) has the following interpretation. 
Let MPN' 1 =1 -P. Then 

v G imP 4$ (I - P)v = <S> PN~ l v = <^> N~ l v G kerP 
v G kerP 44> Pv = O PAf^v = N~ l v O N~ l v G imP 
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Hence, dimimP = dimkerP = d = n/2, which implies dimVi m = dimVk er - 
In this case, given a projector P, one can always find a corresponding J\f. 
Indeed, let v±, . . . ,v d be a basis in imP and Wk := N~ l Vk (k = 1, . . . , d) an 
associated basis in kerP. By virue of Af 2 = a we have N~ l Wk = a^ 1 Vf t . 
Therefore 

A/^Vi,. . .,v d ,W!, ...,w d ) = (wi, . . . , w d , vt/a, . . .,v d /a) 

(where (yi, V2, ■ ■ ■) denotes the matrix with columns v%, V2, ■ ■ •) and, finally, 

Af = (vi, ...,v d ,wi,.. .,w d )(wx, ■ ■ .,w d ,vi/o, . . .,v d /a)~ 1 . (17) 

The TV obtained in this way depends on the choice of the bases V\ , . . . , v d and 
wi, . . . , w d (we can put Avk, det A 7^ 0, in the place of Vk and Bwj, det B 7^ 
0, in the place of Wj). In other words, Af is given up to nondegenerate d x d 
matrices A and B. 

The formulas @ and (|T2*)) coincide after appropriate identification of the 
parameters. Indeed, comparing coefficients by powers of A we have: 



A+ - A_ 

n{-\ + + ( a + - a_)p) = *+<p-*- d -*^ ~_ a-^*^*; 1 

and after straightforward computation we get 

P = (<p + V + d + V+ 1 - (p-ty-d-$>Z 1 )- 1 <p + ^ + d + $>+ 1 , 

I- P = (ip_^_d_^Z 1 - (f+^+d+^fl^cp-^f-d^Z 1 • 
Taking into account the assumption we have: 



(18) 



(19) 



D-D+ -D_D + , N 

D + D_+D_D + cj(A + -A_) 2 ' 1 ; 

The above results are valid for n x n matrix linear problems. Now, we 
focus on the 2x2 case. Because the elements cL are nilpotent (d± = 0), 
then there exist vectors v + , V- such that 

d + v + = , d_v- = . (21) 

Then from (j!9j) it follows immediately Pty+v+ = and (I — P)^_f_ = 0, 
i.e., ^+v+ span kerP and ^f-V- span imP. Hence, v + £ Vter and V- € Vi m . 
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It is not difficult to check that the general form of 2 x 2 matrices d± such 
that d\ = is given by 



d±=[ „T ) = [ T ]( ~a± b ± ) , (22) 



-a±b± b\ \_ I b± 
-a± a±b± J \ a ± / 

where a±,b± are complex numbers. Therefore, to satisfy (|2*T1) . we can take 

—(::)■ ™ 

We have almost unique correspondence (i.e., up to a scalar factor) between 
w + and d+ and between V- and cL. 
Denoting 



we get the explicit formula for P 



B- \ I B+ B_ 



-A + B_ B + B_ 
-A. A_ B+A_ 



P ~ [ A_ J ^ A + A_ ) ~ A_B + -A + B_ (24) 
The corresponding M reads (compare Q17|>): 

_ 1 ( a*-*- ~ A+B+ B\ - aBl \ 

7V " A.B + - A + B„ \ aA 2 _ - A\ A + B + - a AS- ) [Zb) 

Although we can reduce our approach to the explicit formulas (|24[) and 
l|25j) the main advantage of our method consists in expressing the Darboux 
transformation in terms of \I/±d±\I'I|I 1 and avoiding difficulties with parame- 
terizing kernel and image of the projector P which is especially troublesome 
in the Clifford algebras case. 



3 Reductions 

Let us consider the unitary reduction 

t>J(A) = -U>{\) ■ (26) 

If Un is a polynom in A, then the condition (|26|) means that the coefficients 
of this polynom by powers of A are ii(n)-valued. 
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One can easily prove that © implies ^(A)^(A) = C(A), where C(A) 
is a constant matrix (C, v = 0). The matrix C can be fixed by a choice of 
the initial conditions. Usually we confine ourselves to the case 

*t(A)^(A) = k(X)I , (27) 



where k(X) is analytic in A. From (|27j) we can derive k(X) = k(X). By virtue 
of |2J), the Darboux matrix have to satisfy the analogical constraint: 

D\X)D(X) = p(X)I . (28) 

Assuming that D is a polynom with respect to A, compare (JHJ), we get 
that p(X) is a polynom with constant real coefficients, i.e., p(A) = p(X) and 
P,v = 0. 

Lemma 2 If D is linear in A and i28\) holds, then roots of the equation 
det-D(A) = satisfy the quadratic equation p(X) = 0. 

Proof: Let p(X) = aX 2 + (3X + 7. From © it follows 

Aj^o - 7 , - " , A A i + = (29) 

which can be easily reduced to a single equation for S := —AqA^ 1 . Namely, 

aS 2 + 0S + 7 = . (30) 

Therefore, the eigenvalues of S have to satisfy the equation p(X) = 0. Indeed, if 
Sv — uv, then (ap 2 + f3p + 'j)v = 0. On the other hand, the equation detD(A) = 
can be rewritten as 

= det(A/-S , )det^i , (31) 

which means that the roots of det D(X) — coincide with eigenvalues of S. □ 

Lemma 3 ITe assume il(J\) . Then the reduction \2% imposes the following 
constraints on the Darboux matrix fjfi): 

A_ = A+ , dtd+ = , (32) 

and (for n = 2) (v + | vJ) = 0. 
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In particular, by virtue of ©, we can take d- = fd/ + , where / is a scalar 
function. 



Proof: Let us denote zeros of the polynom p(X) by Ai, A2. Because p(X) — p(X) 
there are two possibilities: either A2 = Ai or Ai, A2 are real. From (|10|l we have 

(detD(\)) 2 = a n (X- A+) n (A-A_) n . (33) 

Therefore, in the case (jlOII . Lemma |21 means that A + , A_ coincide with Ai, A2. 

Suppose that A + £ R. Then from (gSJ we have (D(A + )) t Z?(A + ) = which 
implies D + = D(\ + ) = (because for any vector v 6 C™ the scalar product 
(v \ D f + D + v) = 0, hence (D + v \ D + v) = 0, and, finally D+v = 0). Therefore A+ 
(and, similarly, A_) cannot be real. Thus A_ = A+. In this case l|28[) reads 

(Z>(A_))tD(A + ) = . (34) 

Using Q) and lt2*7|) (assuming fc(A±) ^ 0) we get 

(D(A-)) 1 ' = ip-i^flr^l^l = f>-t' + d t _y+ 1 

and llii lfl assumes the form <^ + (p_^' + (iL<i+4'+ 1 = 0. Hence et<i + = 0. 

Finally, in the case n — 2, we use (|22p. Then the condition ct(i+ = is 
equivalent to a + a_ + = 0, i.e., (v+ | «_) = 0. □ 

Another very popular reduction is given by 

[/ M (-A) = JU^X) J- 1 , J 2 = c I , (35) 

then one can prove that Vl/(— A) = J^I / (A)C(A), and we choose such initial 
conditions that C(A) = J , i.e., 

(-A) = J^r(A) J' 1 , D(-A) = JD(\).r l . (36) 

Such choice of C(A) is motivated by a natural requirement that ^, D are 
elements of the same loop group (by the way, the formula (|27|) has the same 
motivation) . 

Lemma 4 We assume U(J\) . Then the reduction Held]) imposes the following 
constraints on the Darboux matrix JIJ); 

A_ = — A + , ip + = ip_ , d + = J~ 1 d-J , (37) 

and (for n = 2) V- = Jv + . 
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Proof: From i|36|l it follows that detD(X) = detD(— A) which means that the 
set of roots of the equation det D(X) = is invariant under the transformation 
A — > —A. Therefore A_ = — A+. Then, using once more we get D- = JD+J -1 
and = Jty+J . Hence <p+d + = ip-J~ 1 d-J. Thus tp+ = Co</?— , where Co 
is a constant. Without loss of the generality we can take Co = 1 (redefining d± 
if necessary). In the case n = 2 the kernels of d± are 1-dimensional. Therefore 
= d + v + = J d~Jv+ implies v_ = ciJv + , where ci = const. We can take 
v + = Jv- . □ 

Other types of reductions (compare |H Ej) can be treated in a similar 
way. 

4 The multi-soliton Darboux matrix 

In this section we generalize the approach of ^3]. First, we relax the as- 
sumption (JSJ. Second, we consider the iV-soliton case (the Darboux matrix 
is a polynom of order N): 



We denote D k = D(X k ), T k = T(X k ), V k = *(A fc ) and U k „ = U^\ k ). 



Evaluating © at A = A^ and multiplying the resulting equation by T k from 
the right we get: 



D(X) = A + A 1 X + ...A N X N . 
The condition © will be replaced by: 



(38) 



D(X k )T(X k ) = 



(39) 



D k ,ti T k + D k U kfl T k = 
To solve the equation (jlHjl we define d k and h k by 



(40) 



D k = y k d k V k 



-l 



T k = $> k h k ^ k 



-i 



(41) 



D k 



i 



Therefore 



and, taking into account (|39|) and (|41j). we rewrite (|4Uj) as follows 



(42) 
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Finally, as a straightforward consequence of ()39|) and (|42j) we get the follow- 
ing constraints on d k and h k : 

d k h k = , d k h k ^ = . (43) 

In ^3] we confined ourselves to the case T(A) = -D(A), i.e, dfc = tp k do k , (ip k 
scalar functions, d$ k constant elements satisfying d^ k = 0), h k = d k . Now we 
are going to obtain the general solution of (|43[) in the case of 2 x 2 matrices. 

Lemma 5 Let d and h are 2x2 matrices depending on x , . . . , x n such that 
dh = 0, dh, fl = and d ^ 0, h ^ 0. T/ien i/iere exisi constants c : , c 2 and 
scalar functions q 1 , g 2 , p 1 , p 2 f depending on x 1 , . . . ,x n ) such that 



d 



12 11 

(j r —qc 

q 2 c 2 -q 2 c 1 

/ / (44) 

11 12 1 

rjr cjr 

c 2 p x c 2 p 2 




Proof: The columns of /i are orthogonal to the rows of d If det(d) ^ 0, then, 
obviously, h = in contrary to our assumptions. Therefore det(d) = which 
means that the rows of d are linearly dependent. Similarly, the columns of h are 
linearly dependent as well. We denote them by p 1 c and p 2 c (where c is a column 
vector). Thus h = cp T , where p T :— (p 1 ,p 2 ). 

dh = means that the columns of h are orthogonal to the rows of d. Therefore 
these rows are of the form q 1 c L , q 2 c ± , where c is a vector orthogonal to c, and, 
finally d = qc 1 - . Thus we obtained l|44|l . 

Taking into account the condition dh,^ — we get 

= qc^ic^p 7 + cp T , tl ) = qc L c^p T =3> c^c,^ = 

This means that c 2 c 1 , M = c 1 c 2 , AI , or c 2 /c 1 is a constant. In other words, c 1 = /c 10 , 
c 2 = fc 20 (/ is a function, c 10 , c 20 are constants. To complete the proof we redefine 
p — > fp, q — > fq, and c k0 — > c fc . □ 

Therefore, 

D(A fc ) = ^(A^gfeC^-^Afc) , (45) 

where are given constant column unit vectors, c k is a row vector orthog- 
onal to c k and are some vector- valued functions (column vectors). We 
keep the notation q k c k = d k , but now in general d\ ^ 0. 

We notice that the freedom concerning the choice of q k corresponds to 
the arbitrariness of the normalization matrix. In particular, the condition 
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imposes strong constraints on J\f . The condition © can be rewritten as 

Qk = VkCk 

The constraint (|39j) implies detD(Afc) = 0. In the case of 2 x 2 matrices 
the equation det-D(A) = (where D is given by (EH1)) has 2N roots (at 
most): Ai, ... A27V. 

Taking any N + 1 pairwise different roots (say Ai, . . . , Aat+i) and using 
Lagrange's interpolation formula for polynomials, we get the generalization 
of the formula ©: 



N+l 

m = £ 

fe=i 



^(A^^c^-^Afc) . (46) 



We have also A r — 1 matrix constraints which result from evaluating the 
formula ifiBT) at \n+2, ■ ■ ■ , ^2N- 

j^g^jr-l^ = ^ ( ^ 

fc=o _ Ao) • • • (Afc — Afc_i)(Afc - A fc+1 ) . . . (Afc — Ajv+i) 

where A = Xn+2, . . . , X 2 n- 
We denote 

Q k := y(\ k )q k , C£ := c^-\\ k ) (48) 

The Darboux matrix is parameterized by 2N constants X k , 2N vector 
functions q k and 2N constant vectors c k subject to the constraints 1)47(1 . 

The crucial point consists in solving the system (|47|) in order to get 
parameterization of the Darboux matrix by a set of independent quantities. 
We plan to express 2N — 2 functions from among Qi, . . . , Q 2 n by other data. 
For instance, we choose Qi,Q 2 as independent functions (they correspond 
to the normalization matrix A/") . 

We rewrite the system (|47|l as 

N+l 

OuvQuC^ + °vkQkCt = , iy = N + 2,...,2N) , (49) 
fc=l 

where 



1 



&vk 



(Afc — A^)(Aa; — Ai) . . . (Afc — Afc-i)(Afc — Afc + i) . . . (X k — Ajv+i) 
1 



(\ u — X\) ... (Ajy — \n){X u - Aat+i) 
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Thus we have a system 1)49(1 linear with respect to Qk- We are going to 
express 2N — 2 vector functions Q3, . . . , Q2N by Qi, Q2 and the other pa- 
rameters: Cfc, Afc. Then, using ((48(1 . we could get 93, . . . ,<?2Ar, etc. However, 
it is better to write 1(46(1 in terms of Qk'. 



N+l 

m = £ 

fc=l 



QkckV-HXk) • (50) 



Taking the scalar product of (fH^l by C\ we get 

Q, = -j: avk f n i\ C '\ Qk, (u = N + 2,...,2N), (51) 

fc=2 <7v0\^u I G l/ 

and the scalar product of vth. equation of 1)49)1 by C M yields 

7V+1 

ErflW = 0> (!/ = iV + 2,...,2iV) . (52) 
fc=i 

This is a system of N — 1 linear equations with respect to Qi, . . . , Qn+i- 
Therefore, we can (for instance) express Q3, . . . , Qn+i in terms of Q\,Q2- 
Then, using ((ST)) , we have Qn+2, ■ ■ ■ , Q2N expressed in the similar way. 

Our method is closely related to the Neugebauer-Meinel approach [H]. 
Let D is given by (|38|). We denote by F(D(X)) the adjugate (or adjoint) 
matrix of -D which is, obviously, a polynom in A. Thus 

D(X)F(D(X)) = w(X)I (53) 

where u?(A) = det(D(A)) is a scalar polynom and I is the identity matrix. 

Therefore, we can put T(A) = F(D(X)) in the formula 1)39)1 and identify 
Afc with zeros of detL>(A). 

In the Neugebauer approach the matrix coefficients A k of the Darboux 
matrix are obtained by solving the following system 

D(X k )y(X k )c k = ,(k = l,...,nN) (54) 

where A& and constant vectors ct are treated as given parameters. Thus one 
has n 2 N scalar equations for (N+ l)n 2 scalar variables. One of the matrices 
Ak, say A^, is considered as undetermined normalization matrix. 
We point out that D(Xk) given by the formula 1)45)1 satisfy ((ST)) . 
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5 The discrete case 



The discrete analogue of Q is the following system of linear difference equa- 
tion 

T^ = U^, = l,...,m) , (55) 

where T u denotes the shift in z^th variable, i.e., (T v ^f){x , . . . , x v , . . . , x m ) := 
\P(z , . . . , x v + 1, . . . , x m ). The Darboux transformation is defined in the 
standard way: 

^> = , T M § = . (56) 

Therefore {T^D){T^) = U^D^, and, finally 

(T M D)[/ M = %D (57) 

If .D 2 (Ai) = 0, then multiplying (|57|) by -D(A) from the right, and evaluating 
the obtained equation at A = Ai we see that the right hand side vanishes 
and we get: 

(T M £>i)l7 M (Ai)£>i = (58) 
where we denote D\ := D{\\). In order to solve (|58[) we put 

where := ^(Ai). Then l|58[l takes the form: 

(^T^xXT^xXT^)^ 1 = . 
Therefore, if 

(T^xH = (59) 

then the equation (|58|) is satisfied. The condition ()59|) can be rewritten (at 
least in the matrix case) as 

Imtii C ker(T At (ii) 

In other words, the sequence of linear operators 

... -» r- 1 ^ -» d x -> -» -> ... 

is an exact sequence [T7] . 

Similarly as in the smooth case we mostly confine ourselves to the sim- 
plest solution of (|59j). i.e., d\ = const which implies d\ = 0. The Darboux 
matrix has the same form Q as in the continuum case. 

Summary. In this paper we developed the approach of ^3] considering 
explicitly the most important reductions, extending our results on the N- 
soliton case, and showing that the discrete case is, as usual, very similar to 
the continuous one. 
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